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In  an  earlier  paper,  it 
was  suqqested  that  peakinq  capacitors 
such  as  are  used  in  nuclear 

electromaqnetic  pulse  simulators  be 
used  to  supply  the  fast  current  rise 
times  which  have  been  indicated  in 
recent  researches  in  the  field  of 
liqhtninq  spherics  and  liqhtninq 
strikes  to  aircraft.  A more 
quantitative  experimental  and 

analytical  study  has  been  undertaken 
to  examine  the  feasilibility  and  the 
siqnificant  parameters  for  such  a 
system.  The  results  of  the  studies 
have  indicated  that  averaqe  liqhtninq 
currents  of  20,000  to  40,000  amperes; 
with  moderately  fast  rise  times  of  100 
nanoseconds  (which  is  an  order  of 
maqnitude  improvement  );  can  be 
achieved  fairly  economically,  but  that 
rise  times  of  the  order  of  30  to  50 
nanoseconds  would  increase  the  cost 
and  desiqn  difficulties  siqnif icantly . 
What  remains  to  be  determined  is  the 
statistical  distribution  of  currents 
and  risetimes  in  strikes  to  aircraft. 
A number  of  fliqht  research  proqrams 
are  currently  underway  which  are 
accumulatinq  data  for  establishment  of 
new  rise  time  test  standards v 

THE  RECENT  WORK  done  in  the 
characterisation  of  liqhtninq  by  the 
Fliqht  Dynamics  Laboratory  and  NASA 
has  revealed  that  the  risetimes  of 
liqhtninq  return  stroke  components  are 
on  the  order  of  100  nanoseconds.  This 
is  an  order  of  maqnitude  faster  than 
the  previously  accepted  risetime  of 
one  microsecond.  The  faster  risetime 
produces  more  enerqy  at  frequencies  of 
aircraft  resonances  than  previously 
anticipated.  Enerqy  couplinq  to  the 
aircraft  interior  can  be  due  to  both 
electric  and  maqnetic  field 
penetration  of  apertures  and  the 
faster  risetimes  increase  this 
couolinq.  At  present,  most  impulse 
qenerators  for  liqhtninq  simulation 
are  desiqned  to  qenerate  one  to  two 
microsecond  risetimes  and  cannot 
simulate  the  newly-revealed  threat. 
Modifyinq  present  qenerators  to 


achieve  the  fast  risetimes  and  large 
voltage  levels  is  not  a simple  task 
because  inherent  generator  resistances 
and  inductances  limit  the  magnitude  and 
risetime  of  generator  output. 


ANALYSIS 

It  is  therefore  desirable  to 
design  a new  type  of  simulator 
suitable  for  testing  aircraft  to  this 
new  threat.  Design  of  such  a simulator 
involves  both  analysis  and  test.  The 
objective  of  the  analysis  phase  as 
discussed  in  this  section  is  to 
provide  analytic  (numerical)  results 
to  support  the  design  of  such  a 
simulator. 

The  analytic  results  fall  into 
three  areas: 

o QUANTIFICATION  of  the  results  for 
a test  cylinder 

o QUANTIFICATION  of  results  for  an 
F-16  in  a proposed  simulator 

o INVESTIGATION  of  pulse  generator^ 
shielding  requirements. 

A test  cylinder  30'  long  and  6' 
in  diameter  is  used  as  both  an 
experimental  and  numerical  test  bed 
for  investigating  the  basic  simulator 
concepts.  Parameter  studies  of  the 
cylinder  response  are  accomplished  for 
variations  in  geometry,  pulser 
wave shaping  elements,  termination 
impedances,  and  spark  gap 

configurations. 

These  cylinder  results  are  then 
used  to  provide  insight  into  the 
response  of  a full  sixe  aircraft  (the 
F-16)  in  the  proposed  full  scale 
simulator  design.  Parameter  studies 
are  also  done  in  this  case  for 
variations  similar  to  that  in  the 
cylinder  case.  It  was  found  that  the 
valves  of  injected  current  and  its 
derivative  were  greatly  dependent  upon 
the  sixe  of  the  peaking  capacitor 
which  was  used.  The  sixe  of  this 


I 


capacitor  is  greatly  determined  by 
economic  factors  but  a modest 
investment  in  a 1000  pf.  capacitance 
can  provide  currents  on  the  order  of 
35-40  KA  with  rise  rates  exceeding 
5X10"  A/m/s. 

One  item  of  concern  is  the 
potential  need  for  shielding  the  Marx 
generator.  This  is  a great  cost  item, 
and  it  is  worthwhile  to  determine  the 
extent  to  which  a shield  is  needed.  It 
was  found  that  construction  of  an 
expensive  quality  overall  shield  is  not 
warranted,  although  an  inexpensive 
isolation  flat  screen  would  be 
desirable.  Signal  to  rise  ratios  of 
45  dB  are  expected. 


TABLE  I 

PARAMETER  VARIATIONS 

o Gap  Locations:  Front,  Back,  Both 

o Characteristic  Impedance:  No.  of 
Hires,  Distance  from  Cylinder,  Size 
of  Wires  (6  Cases) 

c Termination  Load:  Open,  Short, 
Matched , 10  Ohms 

o Scries  Inductance  Ls:  50  nH,  and  for 
30  ns  and  300  ns  Risetimes 

0 Peaking  Capacitance:  1 nf,  e nf,  5 
nf 


RESPONSE  OF  THE  RIGHT  CIRCULAR 
CYLINDER 


MODELLING  APPROACH 


right  circular  cylinder  test 


bed  is  used  to  provide  basic 
configurations.  The  basic  cylinder 
configuration  is  shown  in  Figure  1, 
which  includes  the  Marx  generator, 
measurement  point  locations,  geometry, 
and  spark  gaps  and  terminations.  This 
geometry  is  solved  by  treating  it  as  a 
nonuniform  transmission  line. 

The  model  combines  the  solutions 
for  the  telegrapher's  equations  in  the 
test  fixture  itself  with  the  solutions 
for  the  circuit  which  represents  the 
Marx  generator.  The  solution  is 
accomplished  in  the  time  domain  using 
finite  difference  techniques  (1). 

The  geometrical  inputs  for  the 
model  are  the  per  unit  length 
capacitance  and  inductance  which  are 
simply  related  to  the  characteristic 
impedance  Zo  (2,3).  The  transmission 
line  is  nonuniform  in  that  Zo  is 
different  on  the  end  cones  from  that  on 
the  cylinder.  It  is  not  possible  to 
taper  the  conic  sections  so  that  Zo  is 
preserved,  because  of  the  requirement 
to  have  feed  points  and  terminations 
which  occupy  finite  amounts  of  space. 

The  cell  size  used  was  .0762  m 
and  the  time  step  was  100  ps. 


RESULTS 

A considerable  number  of 
parameters  were  varied  to  quantify  the 
response  of  the  cylinder  in  the  test 
fixture.  The  parameters  considered 
are  listed  in  Table  I. 


RESPONSE  OF  THE  BASIC 

configuration  includes  the  Marx 
generator  of  Fig.  1.,  a peaking 
capacitance  of  1 nf,  a transmission 
line  characteristic  impedance  Zo  of 
71.4  n,  a gap  at  the  front,  a 

termination  impedance  equal  to  the 
characteristic  impedance,  and  a series 
inductance  which  gives  a risetime  of  30 
ns.  The  output  spark  gap  fires  when 
the  voltage  on  the  peaking  capacitor 
reaches  300  kV. 

THE  INJECTED  CURRENT  is  shown  in 
Fig.  2. 

One  may  observe  that  for  this 

case,  the  response  can  be  thought  of 
at  the  sum  of  initial  short  fast 

risetime  pulse  and  a long  slow 

risetime  pulse.  These  two  pulses  are 
relatively  independent  of  each  other 
and  can  be  closely  calculated  by  simple 
analytic  formulas.  Each  pulse  can  be 
thought  of  as  coming  from  an  RLC 
circuit  made  up  of  the  series 
combination  of  the  peaking  capacitor, 
waveshaping  inductor,  and  a resistor 

equal  to  the  characteristic  impedance. 
The  long  pulse  comes  from  the  series 
RLC  circuit  made  up  of  the  Marx 

capacitance,  the  sum  of  the  Marx 
resistance  and  the  characteristic 
impedance.  It  is  easily  seen 
t » 0 


of  the  Marx 
characteristic 
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rigure  3 Simple  RLC  Circuit 
Model 
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that  for  the  simple  overdamped 
circuit  of  Fig.  3,  if  one  defines 

z0  = ATci  \ 


then  the  current  I(t)  is  given  by 


-at  -a. t 

I (t)=V/L(  a+-  a.  ) ) (e  -e  ).(2) 

Thus  each  of  the  pulses  of  Fig.  2 
is  double  exponential  described  by 
equation  2. 

VARIATION  OF  CHARACTERISTIC 
IMPEDANCE 

A parameter  study  of  the  of  the 
characteristic  impedance  of  the 
cylinder  in  the  wire  cage  was  done. 
Wire  spacing  (i.e.,  distance  from  the 
cylinder),  the  number  of  wires,  and 


the 

wire  size 

was  varied. 

The  results 

are 

summari zed 

in  Table  2. 

It  can 

be 

seen 

i that  the 

results  vary 

from  67 

to 

130 

2. 

n,  a span 

of  a nominal 

factor 

of 

1 

Wire 

TABLE  2 

CYLINDER  PARAMETERS 
Wire  No.  of 

Surge 

j 

Size 

Space 

Wires 

Imped 

I 

.1 

5 

16 

71.4 

.1 

5 

8 

102.0 

.1 

5 

32 

67.4 

J 

.1 

10 

16 

109.0 

1 

1 

.1 

15 

16 

130.0 

.6 

5 

16 

64.8 

VARIATIONS  IN  INJECTED  CURRENT 

The  variation  in  injected  current 
as  a function  of  Zo  and  peaking 
capacitance  is  given  in  Fig.  4 and 
Fig.  5 respectively.  It  is  noted  that 
the  rear  gap  case  result ^ in  less 
injected  current,  although  the  exit 


current  would  be  greater  in  this  case. 
Currents  well  over  4000  A are 
possible,  for  a drive  voltage  of300 
kilovolts  and  the  importance  of  having 
as  large  a peaking  capacitance  and  as 
small  a characteristic  impedance  as 
possible  is  evident. 


THREE  DIMENSIONAL  AIRCRAFT 
RESPONSE  IN  A FULL  SCALE  SIMULATOR 


Modelling  Approach 
dimensional  finite 

techniques  (1)  is  used 
the  response  of  a full 


- The  three 
difference 
to  model 

size  aircraft 


in  a full  scale  simulator.  The 
configuration  is  shown  in  Fig.  6.  The 
large  clearances  are  required  to 
provide  sufficient  voltage  stand  off 
such  that  arcing  of  the  aircraft  to 
the  outer  wire  grid  does  not  occur. 
Voltages  exceeding  6 MV  are  expected 
on  the  aircraft. 

The  aircraft  is  an  F-16,  and 
shape  of  the  computer  model  is  shown 
in  Fig.  7.  The  cell  size  is  1 meter 
in  the  longitudinal  direction,  and  is 
1/2  meter  in  the  other  directions.  The 
time  step  is  2 ns.  Because 
approximately  5 cells  are  required  to 
resolve  a wavelength,  the  upper 
frequency  limit  of  the  computation  is 
60  MHz.  The  erected  Marx  voltage  is  MV. 
Measurement  points  include  the 
injected  current  and  E and  H fields 
near  the  cockpit. 


RESULTS 

Rerults  ror  1000  pf  peakinq 
capacitance,  a waveshapinq  inductance 
of  2 uH,  and  a 78  ohm  temination 
resistor  are  illustrated  in  Fiq.  8. 
The  peak  initial  current  of  about  45 
kA  is  unchanqed  with  various 
termination  resistances,  but  the  late 
time  respo  iSes  are  ouite  different,  as 
expected.  The  aircraft  resonates  with 
the  peakinq  capacitor  initially  and 
then  with  the  Marx  qenerator  in  the 
late  time  phases.  The  same  basic 
frequency  of  about  5 MHz  is  aqain 
excited,  as  was  the  case  for  200  pf 
peakinq  capacitance.  Aqain,  the  late 
time  increase  in  currents  for  the  lower 
resistances  are  caused  by  the 
discharqe  of  the  Marx  qenerator.  Other 
oscillations  at  10  MHz  are  evident 
for  hiqher  termination  resistances, 
correspondinq  to  the  half  wave 
resonance  of  the  aircraft  in  the 
fixture. 

Parameter  studies  were  done  for 
various  combinations  of  peakinq 


75-4 


! 


i 


I 

I 


| 

i 


I 

I 


i 


i 


X*-' 


Peak  Current  ( kA) 


capacitance  and  waveshapinq 
inductance.  The  results  are 
summarized  in  Table  3,  which  shows  the 
rise  time  and  injected  current, 
maqnetic  field,  electric  field,  and 
their  derivatives.  The  fields  are 
measured  at  the  at  point  shown  in 
Fiq . 7. 

Several  thinqs  are  evident  from 
this  data: 

The  level  of  peak  injected 
current  is  qreatly  dependent  upon  the 
value  of  the  peakinq  capacitance.  This 
is  primarily  an  economic  decision, 
because  larqe  hiqh  voltaqe  capacitors 
are  expensive.  Because  Ls  probably 
cannot  be  made  much  less  than  8uH, 
peak  currents  of  about  35  kA  are  to  be 
expected  from  a 1000  pf  capacitance. 

VALUES  OF  I OVER  5 X 10"  A/m/ sec 
should  be  obtainable  with  the 
simulator.  This  is  in  the  range  of 
values  preseitly  inferred  from  recent 
LEMP  field  measurements. 

VALUES  OF  E AND  H on  the  same 
order  of  that  caused  bv  nuclear  EMP  (> 
10  V/m/sec  and  10"  A/m/sec 
respectively)  are  possible  to  achieve. 

AIRCRAFT  RESONANCES  are  excited. 
Their  amplitudes  are  of  course 
dependent  upon  the  risetime.  If  the 
minimum  achievable  inductance  of  the 
peakinq  capacitor  output  circuit  is  on 
the  order  of  8-10uH,  then  minimum 
risetimes  ^ n the  order  of  100  ns  are 
probably  the  best  one  can  expect. 


EXPERIMENTAL  STUDIES 

Based  on  the  analytical 
investiqations,  a nominal  system  was 
selected  for  the  low  level 
- experimental  tests  with  voltaqe 

scalinq  of  about  10  percent  and  size 
scalinq  of  approximately  full  scale. 
The  test  voltaqe  was  200  to  300 
kilovolts  usinq  a Marx  aenerator 
i.  drive.  The  resultant  data  was  to  be 

used  alonq  with  the  analytical  data 
to  desiqn  the  full  scale  four  meqavolt 
system.  The  test  arrangement  is  shown 
in  Fiq. 9.  It  consists  of  the  Marx 

. qenerator,  the  peakinq  capacitors,  the 

feedpoint  sparkqap,  all  housed  in  a 
weather  shelter,  the  test  object,  in 
this  case  the  6-foot  diameter  by 

30-foot  lonq  cylinder,  and  the 
downstream  sparkqap  or  termination 
resistor.  Plywood  and  2 x 4 lumber 
were  used  in  fabricatinq  the  test 
arranqement.  The  tests  were  carried 
out  with  the  variation  basically  in 
the  terminations  at  the  far  end  of  the 
j array. 


The  result  of  the  tests  are 
presented  in  Table  4 and  in  the 
oscilloqrams  of  Fiq. 10a,  and  b for  a 
terminated  line  and  a line  with  a far 
end  sparkqap  output.  As  shown  in 
Fiq. 10a,  the  current  waveform  of  the 
experimental  test  arranqement  showed  a 
rise  time  of  approximately  100 
nanoseconds  on  the  linear  portion  of 
the  front  of  wave  followed  by  the  more 
slowly  risinq  pulse  from  the  Marx 
qenerator  drive.  The  time  duration 
of  the  impulse  was  approximately 
three  microseconds.  The  array  spacinq 
was  5 feet  from  the  cylinders  to  the 
wires  and  with  300,300  volts  applied, 
this  resulted  in  a radial  electric 
field  of  about  150,000  volts  per 
meter.  With  the  70  ohm  termination, 
the  current  peak  was  approximately 
4,000  amperes.  With  the  far  end 
output  qap,  the  current  oscilloqram  is 
less  smooth  throuqhout  the  entire 
duration  of  the  wave.  One  of  the 
aims  in  the  desiqn  was  to  investiqate 
the  feasibility  of  usinq  a non-EMP 
peakinq  capacitors  in  order  to  reduce 
the  cost  of  the  system  and  in  the  test 
arranqement,  standard  enerqy  storaqe 
capacitors  were  used.  The  addition  of 
inductance  to  the  input  of  the  vehicle 
was  suqqested  on  the  basis  of  the 
analysis  for  cleaninq  up  the 
wavefront.  The  solution  of  usinq 
additional  inductance  suqqests  that  if 
inductance  is  added  to  the  input  to 
the  vehicle,  then  there  is  no 
requirement  for  providing  special  lew 
inductance  peakinq  capacitors,  and  this 
was  found  to  be  the  case,  but  only  if 
a 100  nanosecond  rise  time  is 
adequate.  If  the  faster  risetimes  of 
30  to  50  nanoseconds  are  required, 

then  the  input  inductance  must  be 
reduced  and  the  capacitors  must  be  the 
hiqher  cost  low  inductance  capacitors, 
at  the  expense  also  of  a less  clean 
wavefront. 


FULL  SCALE  DESIGN 

The  full  scaj.e  test  arranqament 
is  shown  in  Fiq.  6 . The  vehicle 
beinq  tested  is  set  up  on  wood  block 
for  insulation  and  the  array  qrid  is 
supported  by  fiberqlass  lamp  posts.  A 
simple  lumber  structure  is  used  for 
support  of  the  EMP  peakinq  capacitors 
and  these  in  effect  act  as  the  hiqh 
voltaqe  bushinq  to  equalise 
(linearize)  the  voltaqe  equally  across 
the  entrance  to  the  transmission  line 
array.  Hiqh  voltaqe  input  bushinq 
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Flqure  10a.  Oscillograms  for  Terminated  Array. 
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desiqn  is  always  a ditticultproblem 
for  hiqh  voltaqe  systems.  The  peakinq 
capacitors  are  in  effect  also  used  to 
equalize  the  voltaqe  across  the  array 
input  to  help  in  solvinq  this  problem. 

The  output  firinq  qap  and  the 
downstream  qap  or  termination 
resistor  are  fabricated  of  fiberqlass 
qasoline  storaqe  tanks  supported  on  a 
wood  structure.  Simple  plumbinq 
hardware  is  used  for  the  sphere  qaps 
in  order  to  minimize  the  system  cost. 

CONCLUSIONS 

The  analytical  and  experimental 
studies  have  verified  that  the  NEMP 
peakinq  capacitor  approach  can  be  used 
for  liqhtninq  simulators  to  provide  an 
order  of  maqnitude  improvement  in 
current  rise  times  over  existinq 
liqhtninq  simulators  and  with 
relatively  clean  wavefronts.  Also, 
electric  and  maqnetic  field  rise  times 
and  maqnitudes  comparable  to  the 
maqnitudes  measured  in  fliqht  data  can 
be  obtained.  Risetimes  of  the  order 
of  100  nanoseconds  can  be  obtained 
relatively  easily,  but  risetimes  of 
the  order  of  30  nanoseconds  would  be 
more  expensive  and  difficult  to 
obtain.  The  remaininq  Question  will 

be  to  decide  what  rise  times  should  be 
used.  Current  Fliqht  Dynamics 
Laboratory  and  NASA  fliqht  proqrams 
should  help  to  provide  data  for  this 
purpose. 
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